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Abstract 

This letter presents a study of prospects of searching for excited leptons during 
Run II of the Fermilab Tevatron. We concentrate on single and pair production of 
excited electrons in the photonic decay channel in one CDF/D0 detector equivalent 
for 2 fb~^. By the end of Run Ila, the limits should be easily extended beyond those 
set by LEP and HERA for excited leptons with mass above about 190 GeV. 



The Standard Model (SM) of particle physics is known to give results that match the 
current experimental data with a high precision. However, because of well known theoretical 
problems and disadvantages, it is widely believed it can not be a complete theory of elemen- 
tary particles, but rather, likely a kind of effective theory at energies below some scale of the 
order of a TeV. While many models of "new physics" beyond the SM have been suggested 
over the years, one of the most straight-forward ideas proposes that quarks and leptons are 
composite particles. Such a scenario is a recurring theme in nature; molecules are composed 
of atoms, which are in turn composed of nuclei, which are in turn composed of nucleons and 
so on. Furthermore, composite models can explain, in principle, family replication, mixing 
in the quark and lepton sectors, and even make the fermion masses and weak mixing angles 
calculable. 

In most composite models fermions possess an underlying substructure which is char- 
acterized by a scale A [|^] with A about 1 TeV or higher 0. While there is no unified 
model of compositeness, a model-independent effective Lagrangian for excited leptons, orig- 
inally proposed in 0, can be used to model single and pair production in experiments. The 
Lagrangian: 

LiH* = TYi^W. + g'^B^r (1) 

Li*i = ^lcT'^''{f9lw,. + f'9%B,^)h + h.c. (2) 

where / and /' are coupling constants, has been used extensively in a number of phenomeno- 
logical papers which presents ideas on searching for excited fermion production and decay 
to final state gauge bosons in e~^e~, pp, ep, and 67 collisions @, |^, Direct searches for 
lepton compositeness has been done extensively at LEP and HERA each with no 
discovery, but with ever more sensitive limits. Unfortunately, only direct searches for quark 
compositeness have been done at the Fermilab Tevatron [^]. 

In this Letter we present a study for searching for excited lepton production and decay 
for the upcoming Run II of CDF and D0 detectors. We begin with an updated simulation of 
excited lepton production, and optimize the sensitivity by studying the kinematics of excited 
lepton production and decay. With these results we estimate the mass reach and compare 
to the recent results from LEP and HERA. 

To study the mass reach of the Tevatron, we simulated single and pair production and 
decay of excited leptons using the upgraded Fermilab accelerator (1.8 2.0 TeV), and CDF 
and D0 detectors for Run II. The Feynman rules from the effective Lagrangian (Eqn. H) 



are implemented into CompHEP |T0[ using the LanHEP software package. We have 
included into this simulation a complete tree-level calculation which takes into account all 
the spin correlations between excited states production, subsequent decays, and incorporate 
the known NLO corrections. All the partial width and known 2^ 2 cross-section have been 
cross-checked at the symbolic level. Events on parton level generated by means of CompHEP 
have been used as an external process for pythia with the help of the CompHEP-pythia 



interface [12|. The underlying event, jet fragmentation, ISR and FSR are modeled using 



the pythia [ff3l Monte Carlo with the CTEq4l WM structure functions. Since Drell-Yan 



production of excited leptons production is similar to that of Supersymmetric leptons, we 



take the K- factors given in , which only depend on the masses of final particles, and vary 
between 1.23 and 1.24 in the mass range of the search. 

While excited leptons can come in three flavors: e*, /i* and r*, we chose to concentrate 
on the electron since the result of /i* expected to be similar to e*, and r's at the Tevatron are 
still difficult to trigger on and identify. Excited electrons can decay via: e* 67, e* eZ, 
and e* — > Wu channels with the branching ratios shown in Fig. [^. The backgrounds to the 
W and Z decay channels in the hadronic channels are fairly large and the leptonic branching 
fractions needed to identify the W and Z channels are small; however, the backgrounds 
to the photonic final states are relatively small in comparison making them a gold-plated 
signature. 

In the case of the single excited electron production mode, there are two possible signal 
pp e*e e"'"e~7 and pp e*^ — > ez/7. Similarly, pair production gives pp e*~^e*~ 
e"''e~77 and pp — > e*z/* — > 61/77. For simplicity we concentrate on inclusive 667 final state 
for both single and pair production as it covers most of the production in both single and 
pair production cases. Before taking into account the branching ratio to photons, the total 
production cross section is shown in Fig. |^ for the case //A = 10^^ GeV"^ and f = f . 

There are a number of backgrounds to the eej channel. The dominant backgrounds are 
iy7+jets and Z7+jets production. Others include ly+jets and Z+jets, and multijets, where 
jets can fake leptons and/or photons. Studies have shown [jl6l that the fake backgrounds 
can be modeled by the kinematics of the irreducible backgrounds. The backgrounds are 
modeled using the same CompHEP simulation structure described above. Since there is 
an infrared singularity at pj^ = 0, we require Pj' > 10 GeV and ARij > 0.1, where ARij = 



{Arjy + {Aipy, and at the generator level, i and j are any lepton-photon combination. 
We take the K-factor, which has a value of about 1.36 depending on pj^, for the background 
processes from the literature |r7| . 

For both signal and backgrounds, we use a parametric simulation to model the detector 
response. The SHW detector simulation has been shown to be an effective averaging 



between the CDF [jT9[ and D0 ||2^ detectors for Run II. After detector simulation the 



kinematic distributions for both the signal and estimated backgrounds are shown in Figs ^ 
and |[ We find acceptances at about the 0.3 level for signal. 

To maximize our sensitivity we assume, for simplicity, that taking a set of cuts which 
minimizes the expected cross section limit also maximizes our sensitivity to calculate these 
limits. To calculate these limits, we use the signal acceptance and background estimates from 
the simulations above. We add in a single detector with 2 fb~^ of luminosity and we assume 
the experimental parameters of Table |l] for systematic uncertainties. We use a frequentist 
method to incorporate the errors and assume that all errors on acceptance, background 
and luminosity are uncorrelated. To determine the expected limit as a function of a given 
cut using Run II data, under the assumption of no observed signal, the 95% confidence 
level (C.L.) cross section upper limit, for a given cut, is uniquely determined by number of 
events observed in the data. With no signal, the number of events observed in the data, A^, 
fluctuates around the mean number of expected events, M, according to Poisson statistics. 
Thus, an average expected cross section limit is given by: 



Luminosity (fb ^) 


2.0 ± 0.1 


Background error ^{t'"^"'' 


0.1 


Acceptance error — 


0.1 



Table 1: Parameters used to calculate the 95% C.L. cross section upper limits. 



Me* 

(GeV) 


^prod 
(fb) 


B (e* 67) 


Cut 
(GeV) 




expected 
^95% C.L 
(fb) 


expected n 
^95% C.L ■ ^ 
(fb) 


150 


712 


0.400 


140.0 


42.0 


41.3 


16.5 


200 


312 


0.334 


185.0 


18.0 


30.9 


10.3 


250 


160 


0.309 


230.0 


7.64 


23.7 


7.3 


300 


87 


0.297 


275.0 


4.36 


20.0 


5.9 



Table 2: Values for ^ = -r^ GeV ^ optimized points. 



<=T.^9,{N,cut)—-^ (3) 

N=0 

In this way we get an average expected cross section limits as a function of each cut. We 
can then optimize for each mass point as a function of the cuts. 

After a study of multi-dimensional distributions for background and signal, we deter- 
mined that the e^e~7 final state has the best signal to background ratio and kinematical 
distributions. Furthermore, we find that applying a single cut on the kinematical variable 
Mg-y (where e is the electron with highest px) gives the best relation between acceptance 
and N to find minimal value of (T95 for a given mass of the excited electron. This is read- 
ily apparent using Figs. |^ and ^ which shows that the best separation between signal and 
background is the M^^ variable. 

Figure ^ shows the expected 95% C.L. cross section upper limit, (795, as a function of the 
Me-y cut for different masses of e* for f = GeV-^ Placing our cut at the minimum of 
the curve gives our optimization, and final expected cross section limit. These results are 
shown in Fig. Using the Feynman rules for the Lagrangian in Eqn. Q leads to a simple 
relation between cross-section and ^ for signal: 

(Ai)M£,M,0 = (A2)V(£,MeO (4) 
/ _ (T95 

A 10V(10-2,Me*) ^ ^ 



Taking into account this equation, numerical values of signal cross-sections (where ^ = 
10^^ GeV~^) and results for (795 we can make an exclusion plot in the ^ vs. Mg. plane. This 
result is shown in Fig. |^. We compare our results to those of LEP and HERA in Fig. ^ which 
would give the most stringent hmits, to date, for masses above 190 GeV. 



The prospects for searching for excited electrons at the upgraded Fermilab Tevatron are 
excellent. We expect that with a single detector and 2 fb~^ of data should significantly 
extend the mass reach, especially in the low region and large masses. We also point out 
that similar results are obtained for excited muons which significantly improve on the current 
limits which are not producible at HERA. 
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Figure 1: The branching ratios for excited leptons as a function of the mass of the excited 
lepton. 
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Figure 2: Total cross section for the sum of both single and pair production and decay of 
excited electrons in the ee^ final state. 
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Figure 3: A comparison of the kinematic distributions of the final state electrons and photons 
for excited electron and standard model background processes. The figures show the results 
for single and pair production of e*e — > 667 and e*e* — > 667. 
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Figure 4: A comparison of some of the variables studied which gave good separation between 
excited leptons and standard model background processes in the 667 final state. Note that 
the Mg-y variable gives the best separation between signal and background for Mg* = 200 
GeV at Me-y > 185 GeV. 
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Figure 5: The expected 95% C.L. cross section upper limit as a function of the Mg^ cut 
for different masses of e*. There the cross section is defined as a{pp — > ee* — >■ 667) + 
a{pp e*e* — > 667) and £ = 2 fb"^ 
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Figure 6: Optimized expected 95% C.L. for e* eej production using the eej final state 
and 2 fb~^ of data with a single detector. 




Figure 7: Exclusion plot for ^ as a function of e* mass using 2 fb ^ worth of data and 
single detector. 




of data, and one detector, produce the most stringent limits. 
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